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1 Abstract 
O . 

■ We consider quantum effects of a massive antisymmetric tensor field on the dynamics of de Sitter space- 
ft ■ 

^2 ■ time. Our starting point is the most general, stable, linearized Lagrangian arising in nonsymmetric gravita- 
tional theories (NGTs), where part of the antisymmetric field mass is generated by the cosmological term. 
We construct a renormalization group (RG) improved effective action by integrating out one loop vacuum 
• fluctuations of the antisymmetric tensor field and show that, in the limit when the RG scale goes to zero, 

m ! 

£^ , the Hubble parameter - and thus the effective cosmological constant - relaxes rapidly to zero. We thus 
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' conclude that quantum loop effects in de Sitter space can dramatically change the infrared sector of the 
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■ on-shell gravity, making the expansion rate insensitive to the original (bare) cosmological constant. 
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I. INTRODUCTION 



In modern quantum field theories, the vacuum is never really empty. A clear example of this is 
the harmonic oscillator with fundamental frequency uj, whose ground state has an energy E = hw/2. 
This nonzero 'zero point' energy is interpreted as the energy that is always present, even when the 
oscillator is not excited. One could see this as a manifestation of Heisenberg's uncertainty principle: 
the oscillator is never completely at rest l|, |2| . 

The zero point energy of the quantum harmonic oscillator is an example of a very general phe- 
nomenon in quantum (field) theories. The quantum vacuum gets enhanced by loop Feynman 
diagrams. In quantum field theory these diagrams represent a shift in the potential energy and 
such a shift is, like in classical mechanics, in general unobservable. Therefore often the normal 
ordering prescription is used, which removes the zero point energy. However, there are special 
cases, like the Casimir effect, where the vacuum energy does become important. The experimental 



verification 



3!1 of the Casimir effect 



4j however, may not prove the existence of the vacuum energy 



A. The cosmological constant problem 

The problems with vacuum energy start to arise as soon as we start to consider general relativity 
(GR). Whereas in quantum theories, energies are usually observable as the difference between some 
excitation of the Hamiltonian and the ground state, this is not the case in GR. In GR curvature is 
sourced by the energy momentum tensor T^ v according to Einstein's equations 

where R and R^ u are the Ricci scalar and tensor respectively, g^ u the metric tensor, Q is related to 
Newtons constant Gn as Q = (167rGAr) _1 and A 9 is a possible geometrical cosmological constant. 
T^y is sourced by any form of energy, including the vacuum energy and in fact Lorentz invariance 
implies that 

T™ = -(p)g, u , (2) 

where (p) is the energy density associated with the zero point energy. It is thus clear from ([I]) that 
the zero point energy sources gravity in exactly the same way as a positive cosmological constant 
would. In fact we could write 

A c = A p + A g , (3) 



where A c is the total, classical cosmological constant. It is a measurable quantity and current 
observations , [{J, [k]] indicate that 

A c < 3# 2 ft DE ~ 5 x 10~ 84 (GeV) 2 , (4) 

where Hq ~ 1.5 X 10~ 42 GeV denotes the Hubble parameter today and Qbe — 0.73 is the dark energy 
density in the units of the critical density. However, standard calculations in quantum theory tell us 
that 



M 



K ~ -f , (5) 



m 2 



p 



where m p = 1.2 x 10 19 GeV is the Planck mass and Mx is the maximum mass scale which contributes 
to the vacuum energy. In the most optimistic case, this could be the QCD scale (~ 200 Me V), but 
more realistically it may be as high as Mqut ~ 10 16 GeV or even m p . The cosmological constant 
problem can be formulated as the huge discrepancy between the observed value A c and A p . Even 
if Mx is at the QCD scale the difference is 40 orders of magnitude, while if Mx is at the Planck 
scale, the difference is an astonishing 122 orders of magnitude. Of course one could tune A g to 
fit the observations. This however seems very unnatural, in particular when one takes account of 
phase transitions in the early Universe during which the vacuum energy is generally changed by an 
amount typically many orders of magnitude larger than the observed value (j3]). Many attempts have 
been made in the past to try to understand the cosmological constant problem jl|, Q]. Most of these 
attempts depend on new physical ideas, like new symmetries or breaking the equivalence principle. 
However, ordinary quantum field theory has a property that might shed some light on this problem, 
namely the renormalization group (RG) . 



B. Renormalization Group 

It is a well established fact ljj 12| that upon regularizing and renormalizing a divergence 
in quantum field theory, one inevitably introduces a renormalization scale /i. Since this scale is 
arbitrary, physical observables should be independent of this scale. This requirement is imposed by 
the renormalization group equations. Solving these equations leads to the RG improved theory, and 
the coupling constants of the theory, including A, become a function of /i. Since A now runs with 
the energy scale /i, one might hope that at cosmological scales A will run towards zero, independent 
of its huge value in the ultraviolet where (jSj) is calculated. Since the cosmological constant problem 
can be seen as a conflict between scales: the ultraviolet given by Mx where A p is defined and the 
infrared where A e g is measured, it seems mandatory to study the RG behavior of A. 



In [13I [14I [jjj] it is investigated how A runs with /i, when A p is generated by the Standard Model 
particles. A e fj as it is observed is then evaluated at \i — Hq, the Hubble parameter today. The 
result of these calculations is that A varies only very mildly (usually logarithmically) with [i and the 
consensus is that this can never take account of the missing 122 orders of magnitude, though it may 
make the fmetuning issues somewhat less severe. 

Another approach h to study the RG behav 10r of q uantu m gravity Q fl Q. Whi.e this has 
produced very interesting results (discovery of an ultraviolet fixed point), an explicit calculation of 
the effective action is still lacking. The main reason for this lies of course in the fact that we do 
not know what a renormalizable theory of quantum gravity would look like. Another reason is that 
the authors of [17I . [l8| in the course of the running project the effective action onto a certain 



prescribed form (which typically contains A, R and R 2 terms but nothing else). This form may be 
too restrictive since it does not permit e.g. logarithmic dependence on spacetime curvature, which 
plays an essential role in our one loop investigation. 

In this paper we consider in detail the RG behavior of the cosmological constant in a manner that 
substantially differs from the ideas mentioned above. We explicitly calculate the one loop contribu- 
tions to the effective action of an antisymmetric tensor field with a mass given by the cosmological 
constant. Such a field is motivated by nonsymmetric gravitational theories (see section II CI) , but 
our results do not depend on its origins. For simplicity our calculations will be done in de Sitter 
space. Since de Sitter space is essentially empty space, with a cosmological constant, we consider 
our model exemplary for any space where the energy density is dominated by a cosmological constant. 

The Friedmann equation in de Sitter space is given by 

tt9 A r 5An , , 

H 2 = — + — =t , 6 
3 3' V ; 

where 5Aq stands for quantum loop contributions and therefore depends on the renormalization scale 
fi. An important observation of this work is that loop corrections, apart from the \x dependence, can 
dependent on the Hubble parameter (or more generally on the curvature of space-time) in a rather 
complex fashion, 5Aq = <5Aq(/z, H 2 ), such that the self-consistent solution for H 2 in terms of A c as 
specified by Eq. fl6]) is in general of the form, 

# 2 = A cfr (A c ,/i)/3, (7) 



where A c g can be substantially different from A c of the original theory. 

An important question is at what scale \i one should evaluate the quantum corrections. In other 



words: what is our ultraviolet scale and what is the infrared scale where we measure A e //? It is 
often argued in the literature that a natural infrared scale for /i is the Hubble scale. This means 
that the Hubble length H^ 1 is taken to be the longest relevant wavelength. Causality arguments are 
sometimes invoked to justify that. The ultraviolet scale is then given by m p or Mx- We however do 
not agree with this approach. In accelerating space times modes grow 'super-Hubble' and there is 
no reason to assume that these modes do not have an effect on the geometry or are unobservable. 
Therefore a proper treatment of the infrared behavior of the theory must take these modes into 
account and one should let fi run down to the energy of the largest super-hubble mode, which is 
much less then H. As for the ultraviolet scale the problem with choosing m p (or Mx) is that, while 
/i runs to the infrared, the ratio m p / /i becomes huge. This usually spoils the trustability of the one 
loop expansion. Because of these reasons, we define our ultraviolet theory at fi = H, the Hubble 
parameter. This makes sense, since not only is the scale H naturally present in the theory, it is 
also because of the Friedmann equation ([6]) a natural UV scale. This is because the ultraviolet A c 
is huge (see equation ([5])). Furthermore, if we decrease /i, we find that also H decreases. Therefore 
the ratio H/fi stays closer to one then for any other choice of initial /i, which makes the RG results 
more trustable. 



As a last remark, we point out that a substantial difference between A c fj and A c can only be 
obtained if the quantum radiative effects are comparable to the tree level cosmological constant. 
While this necessarily makes the loop expansion problematic, we are not aware of any fundamental 
reason why this compensation should not take place. An example of such a behavior - dubbed 



attractor mechanism - is considered in Ref. 



191 ] . where quantum loop matter fluctuations generate a 



term which exhibits logarithmic dependence on the Hubble parameter, and which compensates the 
tree level cosmological term. The important difference between Ref. 19| and the work at hand is 



that the compensation mechanism in 



19j is realised via the dynamics of a scalar field in an effective 



potential generated by quantum loop corrections, while here the compensation occurs as a result of 
sending the RG scale to zero. 



C. Nonsymmetric gravitational theories 



Nonsymmetric gravitational theories (NGTs) are extensions of g eneral relativity (GR), where the 
standard axiom that the metric is symmetric is dropped [2^] [sj 22 1. One can therefore make a 



decomposition of the general metric c ( 



B 



where (•) and [•] indicate normalized symmetrization and anti-symmetrization, respectively. One of 
the reasons why such a theory is interesting to study, is that it provides a natural, geometric source 
for torsion 23[. Another interesting property is that the massive antisymmetric tensor field might 



act as the dark matter component in the universe [24| 25| 2g|. Since GR is a highly successful it 
makes sense to consider NGT in the limit of a small I?-field. In 27]] it was shown that the only 
Lagrangian that does not lead to an unstable field evolution is given by (for the more precise form 
of this Lagrangian see Eq. (12611 ) 



C 



-g 



Q(R - 2A) + — B 2 - —F 2 
vv ' 144 12 



-m 2 B 2 
4 



0[B l 



(9) 



Here the curvature term R refer to the curvature of the GR background and F^ up is the field strength 
associated with B^. The mass term is given by 



m 2 = wA + —R 



P 



12 144m2 



R 



(10) 



with lj, x an d p undetermined constants from the theory. The cosmological constant, as it enters 
([9]), is the sum of the geometric (A fl ) and matter (A p ) contributions. The fact that the mass of the 
B field is proportional to the cosmological term is extremely interesting, since in principle it induces 
a large back-reaction on the total energy of the vacuum. This is the question we investigate in detail 
in this paper. 

Most of the conclusions of this paper however are not dependent on the fact that we are dealing with 
antisymmetric tensor fields. For example a scalar field with a mass proportional to A, as might arise 
in Kaluza-Klein type theories, would almost certainly lead to similar conclusions. 



D. Overview 



The contents of our paper are as follows. In section [TTJ we discuss some properties of de Sitter 
space and show how to derive the propagator for the 5-field. In section III II we derive the one loop 
effective action and in section IIVI we renormalize the theory and RG improve it. In section |V] we 
calculate the Friedmann equation and show how the £>-field loops alter the effective vacuum energy. 
We discuss and summarise our results in section IVI1 



II. PROPAGATORS IN DE SITTER SPACE 



In this section we first review basic properties of de Sitter space and then we sketch a standard 
derivation for the de Sitter invariant scalar propagator. Having done this, we show how to con- 
struct the propagator for a massive antisymmetric tensor field in de Sitter space. We show that the 
propagator can be written as the sum of two scalar propagators with different amounts of conformal 



coupling. This is similar as in the case of the photon propagator [28|, |29| in de Sitter space. In this 
section we raise and lower indices with Tf v and rj^, respectively, where rj^ = diag(— 1, 1, 1, ..) is the 
Minkowski metric in D space-time dimensions. 



A. de Sitter space 

For pedagogical reasons we begin by introducing four dimensional de Sitter space (later we shall 
work with general D dimensional de Sitter space), which is the hypersurface given by the equa- 
tion [3o| |l9] 

- xl + xl + xl + xi + xl = ±- 2 (11) 

embedded in 5-dimensional Minkowski space-time, where H is the Hubble parameter. The isometry 
group of de Sitter space, 50(1, 4), is manifest in this embedding. We shall use flat coordinates, which 
cover only half of the de Sitter manifold, given by {i = 1, 2, 3) 

1 H 

X = — sinh(iJt) + —xtfe™, 
ti A 



Xi 6 X-i , 

1 H 

X A = — cosh(Ht) - —x^e 11 *, 

— OO < t, Xi < 00. 



(12) 



In these coordinates the metric reads 



= diag(-l,a ,a ,a ) , 



a = e 



Ht 



(13) 



3l| [32] 



Hr] ' 



rj < 



(14) 



which we can write in conformal form by changing coordinates to conformal time 77 defined as adrj 
dt: 

g^ v = diag(-a 2 , a 2 , a 2 , a 2 ) = a 2 ?^ 

We define the de Sitter invariant distance functions 

4 

(15) 



4 1 
Z(X;X') = H 2 VabX a X b > = 1 - -Y(X;X') 



A,B=0 



In conformal coordinates ( |T2l) these functions read 

1 



- ^y(x; x') 



(16) 



z{x\ x') 

y(x; x') = aa' H 2 Ax 2 (x; x') 
Ax 2 (x; x') = —(\f] — v'\ ~ i e ) 2 + I W ~ | 2 ' 

where Y(X;X') = y(x;x'), Z(X;X') = z(x;x'), a = a(r]) and a' = a(r]') are functions given in ( Tl4l) 
and e > refers to the Feynman (time ordered) pole prescription. The function y = y(x; x') is related 
to the invariant length I = £(x; x') between points x and x' as, 



y(x;x 



B. Scalar propagator in de Sitter space 



') =4 sin 2 (~H£(x;x')J . 



;i7) 



The de Sitter invariant scalar propagator for a massive scalar field is the expectation value 



iA(x; x') = (x\ 



-g(D - m 2 - £R 



D 



where Rp is the Z)-dimensional Ricci scalar, which in de Sitter space-time is given by Rn = D(D 
1)H 2 . □ is the d'Alembertian. The propagator (fT8l) satisfies the following Klein-Gordon equation 



y/=g(n -m 2 - £R D )iA(x; x') = i5 D (x - x), 
where 5 D is the D-dimensional Dirac delta. The de Sitter invariant form of ([19]) is 



2 d 2 d m 2 + £R D 



iG(y) 



id (x — x') 
a D H 2 



(19) 



(20) 



where the invariant propagator is defined as iG(y) = iA(x; x') and we used 
d, = (S M z)l = -\Ha(5ly + 2a'HAx,)± 



(21) 



= H 2 (l ~z 2 )-^+ [f^d^z - a~\D - 2)H(d z)) A , 

where z and y are defined in ffTBT) . Eq. (|2"U1) is a hypergeometric equation, whose general solution is 
given in terms of hypergeometric functions 



iG(y) = Cl 2 F 1 



D - 1 



D - 1 



D 1 y 

v\ — ; 1 

' 2 ' 4 



c 2 2F1 



D- 1 



D 



"•f's)- (22) 



where 



(£>-l) 2 m 2 + £R D \ 1/2 
4 i/ 2 



(23) 



The constants c\ and are uniquely fixed if we require that near the lightcone the solution reduces 
to the Hadamard form, while there is no singularity at the antipodal lightcone (which would lead to 



a- vacua 



45J). This means we require: 



H D " 2 ( D 
limiG(y) = . , . „ , n T — 
y^o yy > (4vr) D / 2 V 2 

l\miG(y) = 0, 



(24) 



which uniquely specifies both constants in Eq. (1221) and we arrive at the well known Chernikov-Tagirov 
scalar propagator in de Sitter space |33|] 



iG(y) 



( 4vr )D/2 



D 



D 



— v\ — ; 1 

' 2 ' 



(25) 



C. Calculating the 5-field propagator 



In this section we calculate the propagator for the anti-symmetric tensor field (5-field) in de Sitter 



space. Our starting lagrangian is 
C = 



27] 



( Q(R - 2A) + ^R 2 - ^g^g'PF^F^ - \m 2 g^g^B aP B, v 



(26) 



where R is the Ricci scalar, g = detfg^] and 

m 2 = oo A + xH 2 , Q 
daBp-f + dpB 



1 



(27) 



16ttG 

Notice that in principle there also is a if 4 contribution to m 2 however there is nothing that forbids us 
to set this contribution to zero and we do so for simplicity. We will focus on the B-field contribution 
(remember that we raise/lower indices with rj^/r]^ and work in the conformal coordinates of de 
Sitter space ( fl4|) ) 



12 

-Bypd^ ( a D ~ e 



-B^d"(c 



D-6 



daBf3~ t + dpB 

1 



4 

d-yB a j3 



77 Q VV P - -a D - A m 2 B ap B^ a ^ 



d/j-Bp-y) 



,D-6 



(d a B ia — (D — 6)aHB j0 ^ (d p B^ - (D - 6)aHB°" / 



(28) 



+ 1 -{D - 6)H 2 a D -*B, B°~< - \a D ^m 2 B a pB^rf\ 

where in the second line we have dropped a total derivative. Next we add a gauge fixing term 

1 



46|: 



C 



GF 



--a 
2 



D-6 



(d a B ia — (D — 6)aifB 70 ) (d p B^ - (D - Q)aHB ^ 



(29) 



to obtain 



C B ^ GF = ~ (B^d a (a D ~ Q d a )B M - a D - A m 2 B^B M - 2(D - Q)H 2 a D ~ i B^ B^ 



(30) 



B^V'^B 

2 ^iiv^poi 



where we have defined 



V pa = — 
w 2 



(d a a D - 6 d a - a D -*m 2 )5 [ ;5: ] - 2(D - 6)H 2 a D -%5 [ ffi ] 



(31) 



and 



sit = k - « , = v,j: = v„ + « 



fj. V 



(32) 



The propagator is given by 



iV pa u [ prT A af3 ](x; x') = iria^-q^S (x - x') 



(33) 



We use the following ansatz for the propagator (the subscripts 2 and 3 will become clear shortly) 



per' 



y A a/3 (x;x') = -2a 2 a' 2 fj a [ p fj a ]pA 2 (x;x') + 5f a fj p ][ a Sp\A 3 (x;x') 

With the help of the identity 

(d a a D - 6 d a )(2a 2 a' 2 A) = a D (d a a D ~ 2 d a - 2(D - 3)H 2 )A 

= a D (D -2(D -3)H 2 )A 



(34) 



(35) 



we find that 



(n - 2(D - 3)H 2 - m 2 ^A 2 (x; x') = a~ D i5 D (x - x) 
(n - 3(D - A)H 2 - m 2 ^A 3 (x; x') = a~ D i5 D (x - x'). 



(36) 



So we see that, similarly to the case of the photon propagator in 



that we can write the -B-field 



propagator in terms of scalar propagators with various types of coupling to the Ricci scalar. 



Our A 2 corresponds to Ac of 



281 ] 34j. In fact the propagators A, B and C of these references 



correspond to n = 0, 1, 2 respectively of 

(n-n(D-n- l)H 2 - m^iA^x; x) = a~ D i5 D (x - x). 



(37) 



This equation is solved by fl25]) . with 



(D - 2n - r 



m 2 \ 1/2 



(38) 



iG(y)^iG n (y). 



III. THE EFFECTIVE ACTION 



Since the S-field appears quadratically in our lagrangian (|26|) . it can be integrated out. Thus we 
get (up to an irrelevant normalization constant) 

1 



exp[ir(50] = / VBexpliSig^iB^)} = exp[iS HE (g^)] 



^Det[D 



fivpcr 



and therefore 



T(g^) = S HE {g^) + -Tr ln(P^) + higher loops, 



(39) 



(40) 



where the trace refers to the spacetime integration J d x = j x and the contraction over the Lorentz 
indices. Using f l3Tl) we find 



T B = -Tr \n(D^ upa ) 

= \\ (WM* ^ [ - \ ((d a a D - G d a - a D ^m 2 )rf^ v - 2(D - Q)H 2 a D - 4 5 [ ^ ][p sf ' ' ^ 



The standard technique 12J 35( to get rid of the log is to take a derivative with respect to the 
mass and then write T in terms of the propagator. In our case however, we have to take a log of 
two expressions which are, because of the tensorial structure, orthogonal. To take this properly into 
account we need to take both the derivative with respect to the mass and with respect to H 2 . We 
obtain 



B 



^ 2 (r ?p[ ^ ]ff rA Q/3 ](x;x)(ia^- 4 )^5 ; 



dH 2 ( W , k rA^](x;x)((D - 6)a D - 4 5 [ » V &X 



J dm 2 a D - 1)(D - 2)zA 2 (x; x) — (D — l)zA 3 (:r; a;)) 
+ J dH 2 a D {D - 1)(D - 6)iA 3 (x;x) 
J^a D J dw{D - l)H 2 {^{D - 2)zA 2 (x; x) - iA 3 (x; x) 



-- I a 



+ /^-i,(( D -^„ + (( D -,)-.K^-) 



= Y w + T H , 



(42) 



where 



w 



m 



(43) 



Since at the coincidence limit only the y° term of the propagators contribute (the other powers 
are -D-dependent powers and do not contribute in dimensional regularization) , we use the following 
expression for the scalar propagators. 

t i / n_i \ -n / n — i 



iA n (x;x) 



ttD-2 , 

= JL r(i-- 

/'-tcrm (4vr) D / 2 V 2 



(44) 



with v n defined in (|38|) . In order to renormalize the theory, we need to split the effective action (j42|) . 
with propagators given by (I44I) . in its finite and infinite parts. This separation is done in Appendix IA1 
The final result is 



B 



-9 



\v{\ - D/2) [ 8(D - 5)(D - 6) + 2(24 + D(D 2 - 6D + 4)) X + {D - 3){D + 4) X 2 



4(4tt) d / 2 



2D 



+ H 2 uoA[A + D 2 + D( X - 6) - 



12(X~2) 



+ , ) + ( UA )^- 3 »< D + 4 »| 



647T 



( ((cuA) 2 + 2(1 + x>A# 2 + 2(1 + X + ^)tf 4 ) In ( 



2D 

uA + X H 2 \ 31 
" 15 



A 4 



4 2 v A; 2 v 2 ; 10 



Lff 



SC. 



In this calculation we used the approximation that 

H 2 



<1, 



(45) 
(46) 



cuA + X H 2 
which we justify in section IVB2I 

The effective action (1451) can be simplified further in the special case when \ < 0(1), since then our 
approximation (T46I) is equivalent to H 2 /(ujA) <C 1. In this case we should replace the logs in (j45|) : 

'ojA + xH 2 ^ (uA^ | X H 2 i(xH 2 y 

1]+ ^P_1_(XH 2 ^ (47) 



H 

In 



H 2 
ujA + y# 2 



ln l^ 2 J + c^A 2 V cuA 
/c^Ax y# 2 l/x# 2 y 



IV. RENORMALIZATION 



The Lagrangian we wish to renormalize is (remember that in de Sitter space R = 12H 2 ) 

£ = ^[q{\2H 2 - 2A) + bH 4 + 5£^ 



(48) 



and to do so we add the following counterterms 



C OT = y/=g(A c + Q c H 2 + b c H < 



(49) 



We use a regularization scheme in which the counterterms remove all non-log terms from SC. This 
scheme has the advantage that the equations of motion become relatively simple. The results in 
other regularization schemes are related to this regularization scheme by a finite shift in the coupling 
constants. A potential disadvantage of this scheme is that it is not immediately clear what are the 
physical Newton and cosmological constant. This issue we shall study in more detail in section HVBI 
Using our regularization scheme the lagrangian becomes 



-g[ Q{Y2H 2 - 2A) +bH 4 



64tt 2 

31 
15 



(VA) 2 + 2(1 + X )^H 2 + 2(1 + X + X 2 /m 4 ) In ( 



luA + xH' 



(m 2 ) 



(50) 



Notice that if x ~ ^(1)> Eq. (I46p implies that the replacement (T47j) is applicable, one should only 
replace the logs, since the other terms generated by (1471) are removed in our regularization scheme. 



A. RG improvement 



The renormalized lagrangian (1501) still contains the arbitrary mass-scale \i. The dependence on 
this scale is unphysical, since it is introduced by the counterterms. Removing the n dependence is 



tantamount to the renormalization group improvement 



36| |3j [Jj 



Q 0. 



From 



VdnCrea = (W) 2 + 2(1 + X )^AH 2 + 2 (l + X + ^) & 

we read off the beta-functions 



(51) 



/9a = 


3 (ooA) 2 


1 + x u •A 2 


64tt 2 Q + 


64tt 2 Q 


Pq = 


^-ooA 






64vr 2 










IGvr 2 i 1 + X 





(52) 



The idea is now to 'improve' our lagrangian by imposing the renormalization group equation, which 
in our case (where wave function renormalization can be neglected) can be well approximated by the 
Callan-Symanzik equation 

(fid,i + /3 A d A + PqOq + (3 b d b )C ren = 0. (53) 



We solve this equation by the method of characteristics [37[ [38[ (see also Appendix [Bj), which means 
that we need to make the substitution 



£ rcn (A, Q, b, n) -> C Tcn (A[t], Q[t],b[t],fi[t]), (54) 
where A[t], Q[t] and b[t] are the solutions of the differential equations 

|aM=/5a, j t Q[t]=(3 Q , ^>M=A>, (55) 

and 

H[t] = lie 1 . (56) 
t is a parameter, independent of the couplings A, Q or 6, that we choose to be 

* = ]»(£) (57) 

such that 

fi[t] = if (58) 

This choice of i is however not unique, and different choices lead in principle to different RG improved 
effective actions with different boundary conditions. However, choosing a different t differs from the 
effective action flBTl) only at higher order in the coupling constants, and hence we believe that the 
results presented in this work are generic. On the other hand, since in this work tree level and one 
loop contributions are comparable, this question does deserve further study. 

Our motivation for the choice (j57p is the following: first of all the scale if is naturally present 
in the theory. Moreover, since A in the ultraviolet is huge (see flS}), t = (/i = if) gives a natural 
ultraviolet scale. One could argue that \i = Mx or \i — m p are also natural ultraviolet scales. While 
this in principle is true, such a choice presents the problem that when /i runs to zero, t blows up. 
This is problematic, since this means that the running parameters hit the Landau pole (see fl59|) ). 
Our choice (15 7p on the other hand does not present this problem, since we will show that if decreases 
when /i decreases in such a way that t never blows up. Thus the choice (15 7p has the advantage that 
it results in trustable running from the far ultraviolet to the far infrared. 
The solutions to the renormalization group equations (1551 are given by 

i+x 

Q[t] = Q (l - 2(2(1 + x ) + 3u;)Qo WJ 

A[t] =A„(l- 2(2(1 + x )+3w)Qi"W) < 59 ) 

w= 6 °-iip( i+ * + k)<< 



where a subscript zero means evaluation of the parameter at t = 0, and 



Q 



= 1287r 2 Qo = 8vrmJ, 



(60) 



where m p = 1.2 x 10 19 GeV is the Planck scale. 

With the replacement (I54|) we get the final RG improved lagrangian 



C 



RG 



64vr 2 



Q[t](R - 2A[t}) + b[t}H 4 



uA[t}) 2 + 2(l+ X )coA[t}H 2 



(61) 



2 (-^ + ^^)^ 4 ) ln ( 



H 2 



This effective lagrangian is finite and independent on the scale fi (in the sense of the Callan-Symanzik 
equation), and we use it in the following to study how the one- loop B field quantum corrections 
influence the Friedmann equation in de Sitter space. 



B. Physical parameters and boundary conditions 

Since H appears in the logarithm of flBTl) . it is a priori not clear what we mean with our physical 
parameters. Since we - naturally - do not want our model to spoil e.g. solar system measurements, it 
is important to identify physical parameters. We define the physical constants (denoted by subscript 
p) in a standard way 

1 C 

A[t] p Q[t] p = -- lim 



2 H^O y/=g 

Our analysis in section IVl shows that A p [t] is not the relevant parameter that determines the expansion 
rate of the universe. Instead, the relevant quantity which determines the expansion rate is given by 
A off , which is the self consistent solution to the quantum Friedmann equation. 

In the following we assume uj x + 1 an d x < C(l)> such that we can apply (147|) (see also the 
note after equation (150]) ). In this limit, the RG improved parameters fl59|) become 

Q[t] ~ Q 

A[t] ~ l-e^fOd-iAnt (63) 



6uj 2 (Qq) 



With these parameters, we calculate ( 162]) for the lagrangian ( IBTj) and obtain 



Q P [t]A P [t] 



lim 



Q[t]A[t] 



3 2 »r l2n fuA[t] 

— ^ 2 A[t] 2 ln(-^i 



1287T 



2Q o A[t]+0(^,tf 2 



(64) 



To calculate the first derivative of Crg/ y/—g, first note that the derivative of the parameters Q[t], 
A[t] and b[t] can be easily obtained by noting, dt = —f^d^, which is evaluated in Eq. ( 15T1) . whereby 
the potentially divergent 1/H 2 terms cancel out. The result is, 



Q P [t] = Q[t] 



Qo 



where we used 



uA[t] 

1287T 2 

_ 1 + X 
3a; 



;i + x)in 



a; Aft] 



(65) 



2t 



-2t. 



(66) 



Equation fl65|l implies that Q p is equal to Qq plus a small, t- independent, correction (remember that 
we assumed ui ^> x + !)• Since this means that the Newton constant does not run with the scale fi, 
standard gravitational tests are not affected in this limit. 
Furthermore we find that, since Q p — Qo, Eq. (|64j) implies that 



A p [t] ~ 2A[t] 



(67) 



and therefore it makes sense to put (see Eq. (jSJ)) 



A = M x A /m 2 



Unfortunately it does not seem to be possible to make a similar statement for b. However this is of 
no great concern to us, since the term b^H 4 does not contribute to the Friedmann equation and thus 
it is not expected to influence gravitational tests. 



V. THE FRIEDMANN EQUATION 

From 059p it is clear that A runs only logarithmically with \i and this is never sufficient to get 
A[/i] — > in a satisfying manner Q]. The physically relevant quantity however is the effective 
cosmological constant given by the (modified) Friedmann equation 



where A c fj may be very different from A. To get the Friedmann equation we need to calculate 
<9£ RG 



dg 



° ^ "2 £rg/ ^ + 4^ dtf 2 " + 8 dt 







(70) 



for the RG improved lagrangian fIBTl) . After some algebra and dividing the result by — 3Q we obtain 



i^i + (i + x)^l 

A[t] f 3u 2 A[t] 
1 



Q[t] J Q[t] 
1-(1 + X + 3^ A ^ 3 ^ 



;i+x) ln( 



uA[t] + 



H 4 
+ -~r. 



Q[t] 



1 + X+ 2 X+ 3 



(1 + X + 3a;) 



Q[t] J Q[t] 

u 2 A[t] 



In 



# 2 

^A[t] + X g 2 
# 2 



(71) 



where we have removed a H A /m 2 term, in agreement with our approximation (|46[) . In order to find 
the effective cosmological constant, we need to find the self consistent solution of (171j) for H 2 . 



A. Case 1: uj 2 K[t\/Q[t] <C 1 

At first instance, this case might appear to be the most interesting to look at. In this limit the 
Friedmann equation reads 

tt2 Alt] H 4 31 1 9 

H 2 ^ + a^- = a = — + x + ~X , 72 

3 Q[t] 30 A 2 A K J 

which is solved by 



H 2 = Q[t]^ ^ . (73) 



Clearly there is no way in which (J73j) would lead to A cff — >• 0. To compare (JT3j) with the standard 
form of the Friedmann equation, we assume aA[t] <C Q[t] and get 



H 2 = m 



aA[t] , G ^A[t] A 2 



J 



3Q[t] \3Q[t] 

Therefore we see that we get a small correction to the Friedmann equation. 

On the other hand, one could assume that aA[t] ^> Q[t}. In this case Eq. ( 1731) implies that H 2 



(74) 



Q[t]A[t]/(3a), such that H 2 is given by the geometric mean between A[t] and the Planck scale 
Qo = 87rm 2 . This clearly does not help us in addressing the cosmological constant problem. Moreover, 
since in this limit perturbative calculations are expected to break down, this result is not to be taken 
too seriously. 



B. Case 2: ujA/Q < 1 and uj 2 A/Q ~ 1 



This limit requires that Mj is considerably smaller then m p and u needs at least to be order 10. 
However these assumptions are not unrealistic. For simplicity we also assume that u % + 1, so 
we can use the results of section HVBI Furthermore we assume that \\\ < 0(1), so we can use the 
approximation (|47[) (see also the comment following equation (]50p ). Finally we shall drop the H 4 /Q 
terms in the Friedmann equation, since we will be interested in the regime where H becomes small. 
We now summarize all of our assumptions, 



H A /Q < H 2 < uA < Q , w 2 A~Q, cu>x + l |x| < 1 • 



(75) 



With these assumptions the Friedmann equation (1711) becomes 



H 2 



where we have defined 



A[t]Q 

3u 2 
X 



9 



3A)3AL- |A 2 



(76) 



u 2 A[t] 
Qo 



6Xnt 



u 2 A 



o u 2 M 4 x 



57T mz 



(77) 



In 



QoA 
uoH 2 



In the second line of ( 1771) we have used the results of section IIVB1 

If we start the RG flow at \x = H and let \x go to zero (see section [IB]) , there are two possibilities. If 
H decreases slower then /x, X[t] increases; otherwise X[t] decreases. Furthermore it is clear that if X 
increases, and H decreases, clearly L will increase. We solve the equation [r.h.s. of (I76p ]=0 in this 
limit and get (apart from the trivial solution) 



> 1 11 11 ( 1 



(7f 



Therefore we see that when H goes to zero, A approaches 1/3 from above. But this means that A 



ought to decrease. Fortunately this is exactly what happens: for A < A crit |47|, the ratio H/fi grows 
and thus A increases. However when A equals A cr it, H starts to decrease faster then /x, so from that 
point onwards, A actually starts to increase, driving H 2 towards zero. Before studying this process 
in more detail in section IV B 11 let us first determine how fast H 2 goes to zero as A \ 1/3. To do 
this, we rewrite the Friedmann equations (176|) as follows: 



H 2 



QoA 



10 



exp 



3A)3A 



9 



A^ 



^ 2 H 
QoA 



2s 1 



(79) 




The limit H — > and A \ 1/3 from above thus means approximately: 



tt2 Qo 

H ~ — exp 

3a; 



(80) 



6(A- 1/3) _ 

From this equation we see that, since A depends on \n(H/fi), we actually have approximately that 
H oc fi and thus we get power-law running instead of logarithmic running. In fact H goes even faster 
to zero, since we require that H goes to zero faster then \i and therefore A 'accelerates' towards 1/3. 



1 . Numerical analysis of the model in the infrared limit \i — > 

In this section we demonstrate that the behavior described in the previous section, is actually 
realized in this model. Unfortunately we are unable to solve the Friedmann equation analytically 
for H, therefore we need to rely on numerical analysis. For definiteness we choose Qq = 1000 and 
uj = 10, however none of the qualitative features of the model depends on these numbers. One does 
not need to specify Ao, since Ao will only define the 'starting point' on the curve ( 1VB II) . While we 
choose two parameters, we actually have a one parameter family of curves, since the Qq dependence 
in f|75|) can be removed by defining h 2 = H 2 /Qq. 

We solve the Friedmann equation numerically and plot H 2 versus Aft] (figures IV B 1 1 and |2|) . The first 
thing we need to consider is whether, if we start with < Aq < 1/3 and let /i decrease, \[t] increases 




FIG. 2: Closeup of the region X[t] > 1/3. The solid curve is the solution to the Friedman equation. The 
dashed curve represents Eq. (f80]) . 



or decreases. To see what happens, we consider the Friedmann equation in the limit X[t] — > 



X[t]Q 



3a; 2 

and take the derivative with respect to fi of this equation to obtain 
OH 2 / Q X[t] 2 \^2Q X[t} 2 2Q X[t} 2 



H 2 ~ (81) 



, 0{X[t]f < (82) 

Since this expression is always negative, we find that upon decreasing /x, H will increase and therefore 
X[t] will increase. In other words: starting at the origin of figure IV B 11 we run along the curve to 
the right. It is important to realize that we cannot change the running 'direction' on the curve. Also 
'stopping' on the curve is not possible as long as we keep decreasing \i. 

Depending on the value of A (and thus Ao)we see that one actually starts with an effective cos- 
mological constant that is somewhat larger then one would expect from the unmodified Friedmann 
equation (shown by the dashed line in figure IVB lj) . Continuing along the curve, at a certain point 
H starts to decrease. In fact, it is clear that H decreases faster and faster until a point where it 
decreases faster then \i and therefore eventually Aft] starts to decrease. This must be the case, since 
this is the only solution to the Friedmann equation in this regime and we see no good reason why the 
Friedmann equation should break down. In figure [2] we show how H goes to zero, compared to our 
estimate flHUj) . We see that H decreases even faster then exponentially as X[t] approaches A cri t = 1/3. 



The closer X[t] to 1/3 is, the more accurate (IHUi) becomes. 

Notice that the fact that A stays finite for all values of /i between H and is an important property 
of our t parameter (|57|) . For other -at first instance sensible- choices of t like t = \n(M x //j,), A and 
therefore A will hit the Landau pole become infinite. 

2. Does H 2 stay smaller then ujA ? 

In calculating the effective potential we assumed that (146]) holds, or in the present case, where 
Ixl < 1: 

H 2 < uA. (83) 

We shall now show that this assumption is indeed justified over the whole range of X[t], independent 
of the parameters. 

First of all notice that the assumptions f!75|) mean that 



u > 0(10). (84) 

Therefore it is clear that if the standard Friedmann equation, H 2 = A/3, would be correct, the 
assumption ([8*31 holds. To see whether this is also true for our modified Friedmann equation, we take 
the derivative of ( |76l) with respect to A 

w ^ (3A - ^ + i ( 2 + 6A - 45A2 + 6 < 2 - 9A > Ai ) • (85 > 

with L defined in (177j) . To calculate at which A H is maximal, we put L approximately constant and 
find 

dH 2 2L + 1 + V4L 2 + 16L + 11 

^T = A = s^TT^ (86) 

and we substitute this in our Friedmann equation, divided by ujA to obtain 

H 2 (L + 1) \2L (2L + v/4L(L + 4) + 11 + 12) - y/AL(L + 4) + 11 + 19 

a;A 



7) 



3(6L + 5) 2 cu 

Since we know that in this regime H 2 = O(A), it is safe to say that L grows to as most O(10). With 
this estimate and ( 1841) we get approximately 

IP 

ujA 

and thus our assumption in is satisfied. 



~ 0.1 

max 



3. The infrared sector 



If the mechanism described above is indeed responsible for the small effective cosmological constant 
measured today, we require that \x runs at least to a value /i m j n such that H(n m i n ) = Hq. Here 
Hq ~ 1.5 X 10 -42 GeV denotes the Hubble parameter as it is measured today. The analysis presented 
in this subsection is oversimplified, since we are neglecting matter and radiation contributions in 
assuming a pure de Sitter space-time. We can estimate /i m i n by inverting the relation (150]) , 



H exp 
Hq exp 
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(89) 



3cj 2 A 

Notice that the subscript zero for H means 'today', while for A it means 'at fi = 0. 

According to the logic in section (II Bp the scale (189]) corresponds with the largest scale k amplified 

during primordial cosmic inflation. Since the number of e-folds is given by 



tot 



(90) 



with Hi the Hubble parameter during inflation, we get 



tot 



60 + In 



60 



in nip 
9u 2 H 2 



(91) 



where we used the usual estimate that In — 60 and we assumed that Hj ~ Ao/3. We can 

rewrite (191]) by noting that the power spectrum of curvature perturbation has an amplitude squared, 

H 2 

-2.4xl0" 9 , (92) 



A 2 ~ 



ixera 
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where e is the slow roll parameter in inflation, defined as e = —(dHj/dt)/H 2 (in scalar inflationary 
models its numerical value is typically of the order 10 -2 ). With this the total number of e-foldings (]9T]) 

becomes, 



tot 



60 + 
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60 + 2 x 10 8 x 



uj 2 e 



(93) 



This provides an estimate of the duration of inflation in our model, which decreases as 1/u; 2 as u 
increases. Even though we are unable to make a direct estimate of the residual dark energy density 
in our model, in Eq. (193]) we have provided a link between the the dark energy density (the Hubble 
parameter today) and the duration of inflation. 



VI. DISCUSSION AND CONCLUSION 



In this paper we consider the one loop vacuum fluctuations in de Sitter space of a nonsymmetric 
tensor field, with a mass generated by a coupling to the cosmological constant and the Ricci scalar. 
While such a field is motivated by nonsymmetric gravitational theories, the present discussion is 
independent of this fact. We believe that qualitatively the results presented here also hold for e.g. a 
scalar field whose mass is generated by A. The calculation presented in this work shows that, as one 
would expect, the coupling constants of the theory A, Q = (IGttGat)" 1 and b - which correspond to 
the coupling constants of the terms R°, R and R 2 in the gravitational action, respectively - become 
renormalized. The renormalization depends, of course, on the renormalization scale fi and after RG 
improvement we obtain the standard result that the constants (in principle) run logarithmically with 
\i. Usually this logarithmic running is used to argue that the RG flow cannot drive the cosmological 
constant to zero. Indeed it seems implausible that such a mild running could explain the smallness 
of AGn ~ 10~ 122 , which defines the gravitational hierachy problem. 

The main thrust of the present work is the observation that the observed, effective A g may dramat- 
ically differ from the (RG improved) parameter A. The reason is that A c g- is given by the Friedmann 
equation, H 2 = A e g/3, and this equation gets modified by the loop effects. This modification is 
typically of the order m 2 /m 2 , where m is the mass of the field whose contribution we are calculating. 
This quantity will therefore be in general suppressed. However in our case m ~ A and standard 
quantum field theory calculations tell us that A may be comparable to m 2 . In other words, a large 
A implies substantial modification of the Friedmann equation by loop quantum corrections. At the 
same time our analysis shows that the Newton constant remains to a good approximation frozen. 
This is important, since a fast running of the Newton constant could be disastrous in the view of 
Solar system tests of GR. 

Furthermore, we find that as fi runs towards zero, H - and thus A e fj - changes mildly for some 
range of \x. However from a certain value of /i, H starts to decrase rapidly, eventually even faster 
then linearly, H oc fi. Therefore, while A changes only logarithmically with /i, A e g- can decrease as 
power law. In this manner our model circumvents the usual statement that RG flow cannot relax 
the cosmological constant. It is important to stress that this is independent of the size of the initial 
A. The initial conditions influence only at what value of /i the dramatic decrease in A e fj begins to 
take place. 

An important point we wish to address now is why do we sent fi to zero? To understand this, 
it is first of all necessary to understand that the starting point of our running, \i — H, is in the 



ultraviolet, where both A and H are of the order the Planck scale m p . However in de Sitter space the 
Hubble parameter is not the lowest possible scale. In fact, many physical modes have wavelengths, 
that during an inflationary de Sitter phase get amplified and their wavelength grows larger then the 
hubble length, H~ l . In order to get a theory that takes proper account of these super-Hubble modes, 
the infrared \i should be taken as small as the smallest energy of these modes. Since we are interested 
in the behavior of the theory on large (cosmological) scales, this is exactly what one should do. 
Another interesting question is how applicable our model is to the real world. Strictly speaking, our 
results are only valid in de Sitter space and one could wonder if our results also apply to, for example, 
quasi-de Sitter spaces where if is a (mild) function of time. Quasi-de Sitter spaces are relevant for 
inflationary models and also in the case whenever A e g dominates the energy density of the Universe, 
as it is today. Even though we have performed our analysis in de Sitter space, we expect that our 
analysis should apply whenever dH/dt <C H 2 . This condition is satisfied if A e ff/(87rGAr) dominates 
the energy density of the Universe, which is also (marginally) true today. 

There might be implications of our model for inflation. One could envisage a scenario, where 
inflation is caused by a large cosmological constant and as times goes on more and more of infrared 
modes are produced, pushing the relevant RG scale \x further and further to the infrared, such that A c g 
eventually decreases, terminating inflation, resulting in a standard decelerating Friedmann universe. 



This type of reasoning has also been pursued in Refs. 39 



40 



41] by using different techniques. An 



explicit calculation is necessary in the present or similar model however in order to illuminate the 
questions concerning inflationary dynamics and cosmological perturbations. Based on the assumption 
that the scale \x corresponds to the largest scale amplified during inflation, in section IV Bl we derive 
a relation between the total number of e-foldings during inflation in our model (driven by a bare 
cosmological constant) and the dark energy density as measured today. However at the moment we 
are unable to provide a direct estimate of the dark energy density. 

While our model shows interesting behavior, there are some potential problems. First of all one might 
argue that our model suggests that at small scales A e g is huge, and this could present problems to e.g. 
solar system measurements. However, this reasoning is not correct. Even in small scale experiments, 
the large scale modes still have their effect on the vacuum. In other words: it is impossible to decouple 
the infrared modes from the vacuum. 

a more serious concern is the validity of the perturbative approach. Indeed in our model tree level 
and one loop contributions are comparable and therefore it is not clear whether our results hold if 
one includes higher loop corrections. Furthermore there is a similar problem in the renormalization 
procedure. The choice of the renormalization scheme and the running parameter t are in principle 



arbitrary. Different choices however differ only at higher loop order. Since higher loop contributions 
are probably substantial in our model, there is no way of telling which effects are induced by our 
choices and which effects are really physical. 

Despite these problems we do believe that the following is true. A field whose mass is generated by 
the cosmological constant might generate a huge back-reaction on the vacuum. This back-reaction 
can be of the same order of magnitude as the original cosmological constant. It might be that the 
exact structure of the modified Friedmann equation is such that, because of the RG running of A, 
the effective cosmological constant becomes tiny, without fine tuning and regardless of the initial 
condition. In this work we have shown an explicit example of this behavior. 
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APPENDIX A: SEPARATING THE FINITE AND INFINITE PARTS OF Tb 



To separate the infinite from the finite part of the effective action T (j42p . we use the following 
expression 
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where v n = v n \ D=i and ip(z) is the digamma function, defined by 
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1. Evaluating Y w 

By making use of Eqs. ( 1421) and ( 1A1I) we can break T w into the infinite and finite parts as follows: 
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This breakdown is not unique (since we can always shift a finite part of the infinite part back 
and forth), but when taken together with the counterterms, this nonuniqueness has no physical 



consequence. The integral in the infinite part is easy to evaluate: 
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To evaluate the finite part we assume that 
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which we justify in section IVB2I With this assumption we can use 
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2. Evaluating T# 



From Eq. (|42H we can write 
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We start with the infinite part: 
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The finite part is easily evaluated and gives 
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3. Combining the results 

Since in D = 4 our counterterms are of the form R°, R 1 , R 2 ..., we need to expand the D dependent 
powers in H . To achieve this, one needs to introduce an arbitrary mass scale /i, as can be seen from 
the following identities: 
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Notice that this statement is equivalent to the argument used in e.g. [121 ]. where /x is introduced in 
order to give physical constants the right dimensionality. Using (1A10I) . we combine Eqs. (IA3I) . (I AT)) . 



( 1A8I) and ( 1A9I) to obtain our final result for the 1-loop effective action for the S-field 
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+ X +(c^A) : 



: (D-3)(D + 4)V , 



2D 



+ 



64tt : 



((M) 2 + 2(1 + xVAtf 2 + 2(1 + X + ^)^ 4 ) In (^#^) " ^ ( 



ujA + xH' 
IP 



4 2 v A; 2 v 2 ; 10 



tt2 



= / <5£ 



(All) 



APPENDIX B: METHOD OF CHARACTERISTICS 

In this Appendix we shall describe how the method of characteristics works [37[[38[. The require- 
ment that physics is independent of /i, is equivalent to requiring that the theory is invariant under 
the transformation 

/i — > pte* = p,[t], (Bl) 

with t a real number. Since our coupling parameters A, Q and b depend on /i, the transformation 
(IB II) means that 

A -> A[i] 

Q Q[t] (B2) 
&-►&[*]. 

Quantities that we require to be independent of fi, like the effective action, need to be independent 
of this transformation and therefore 



T(A,Q,fi) = T(A[t},Q[t],b[t},m)- 



(B3) 



This can only be true if 



Comparing this with the Callan-Symanzik equation 

K +p ^ +p ^ + ^) r(A ' Q ' b ^ )=o ' <B5) 

we conclude that the transformation ( 1B1I) and ( 1B2[) . with parameters given by 

dA[t] 



dt 
dQ[t] 
dt 

db[t] 



dt 

enforces the Callan-Symanzik equation for the effective action F. 



/3 A [A[t],Q[t\] A[0] = A 
{3 Q [A[t}} Q[0]=Q (B6) 

A, 6[0] = 6 
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